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^-|s , Precision spectroscopy on cold molecules can potentially enable novel tests of fundamental laws 

fvj ■ of physics and alternative determination of some fundamental constants. Realizing this potential 

requires a thorough understanding of the systematic effects that shift the energy levels of molecules. 

We have performed a complete ab mitio calculation of the magnetic field effects for a particular 
^ . system, the heteronuclear molecular hydrogen ion HD"*". Difi'erent spectroscopic schemes have been 

■^^ ' considered, and several transitions, all accessible by modern radiation sources and exhibiting well 

controllable or negligible Zeeman shift, have been found to exist. Thus, HD"*" is a candidate for the 

determination of the ratio of electron-to- nuclear reduced mass, and for tests of its time-independence. 
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Oj, Cold molecules have recently been proposed as novel systems for precision measurements related to fundamental 

aspects of physics, such as the measurement of the electron-to-nuclear mass ratio, its possible time-variability, tests 

of Lorentz Invariance, tests of QED, and parity violation. Several of these are based on precision spectroscopy [l|- 

O • ISJ, where frequencies of ro- vibrational transitions exhibiting high quality factors must be measured. As the tests 

jrt ' mentioned above have already been performed with very high precision using various atomic systems, studies on 

molecular systems must be conceived in ways that have the potential of surpassing atomic tests. This implies that 

^ ' molecular systems must be experimentally accessible and that systematic shifts of the transition frequencies must be 

' ^ [ sufficiently small so as to guarantee the desired spectroscopic accuracy. 

^^' We focus on a particular diatomic molecule, HD"*". Because of its relative simplicity it can be analyzed with high- 

f~i '■ precision ab initio methods. QED calculations of ro- vibrational frequencies have reached a relative accuracy of a few 

Qh' parts in 10""'^° Q, and the influence of external perturbing and exciting electromagnetic fields, i.e. all systematics, 

can also be treated accurately ab initio. Comparison of theoretical and experimental transition frequencies in HD+ 

fvj , can potentially lead to determining the electron to proton/deuteron mass ratio and the Rydberg constant in an 

^ ' alternative way. For these purposes, an experimental transition frequency accuracy in the range 1 • 10~^° to 1 • 10~^^ 

0^ , must be attained. The requirement is even more stringent, at the 10^^^ level, if the system is to be used to test the 

Cn ■ time independence of the mass ratio or Lorentz Invariance. Experimentally, HD"*" can be cooled to tens of mK by 

^^ ' sympathetic cooling in an ion trap [l^l and rotationally cooled [ll[ . One-photon laser spectroscopy of ro- vibrational 

; ,. transitions has been performed and one transition frequency has been determined with a relative uncertainty at the 

t^^ ■ 2 • 10~^ level Q, in agreement with the QED calculation [9]. 

^^ , Here we present the results of a thorough study of the magnetic field effects on the radiofrequency, rotational 

(THz) and ro- vibrational transitions in HD+. Both one- and two-photon transitions are of interest for precision 
spectroscopy. Earlier, Karr et al. [12l| had evaluated the two-photon transitions strengths in HD+ in the spinless 
particle approximation. An analysis of a particular ro- vibrational transition at moderate spectral resolution has been 
reported in Ref.Q. 

Theoretical approach. The Hamiltonian H of the HD+ ion in an external magnetic field has the form H — H^^ + 
^ . ydiag _|_ ^spm _|_ ^mag ^ whcrc H^^ is thc non-relativistic 3-body Hamiltonian and the correction terms collect the spin- 
independent interactions, the spin interactions (cf. Ref. |13l. Il4|) and the external magnetic field interaction terms, 
respectively. In the leading order approximation T/^'^s = — B • ^j(eZi/2Mic)(Li -I- {fii/si)Si), where the summation is 
over the constituents of HD+ (i = p, d, e), Mi, Zi and /i^ are the mass, the electric charge (in units e) and the magnetic 
moment (in units eh/2Mic) of particle "i", Sj is the size of its spin (1/2 for i = p,e and 1 for i = d), L.^ and S^ are 
the orbital and spin angular momentum operators in the center-of-mass-frame, and B is the external magnetic field. 
Neglecting the higher order corrections to V"^'^^ (cf. Ref. [l^) is justified for magnetic fields B below the threshold 
Bthr ^ 10^ G for which the contribution of V™'^^, increasing with B, reaches the order of magnitude of the hyperfine 
energy, since the relativistic corrections to V"^^^ are smaller than the theoretical uncertainty of the hyperfine energy 
levels of Ref. [3. 

The spin structure of the ro- vibrational state (w, L), where v and L are the vibrational and total orbital momentum 
quantum numbers, is calculated in first order of perturbation theory using an effective Hamiltonian V that is 
obtained by averaging T/^^p'" + \/mag ^^^j. -j-j^g spatial degrees of freedom. Compared to iJcff (the HD+ effective spin 
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Hamiltonian of Ref. [11]), V^ includes 4 additional terms, originating from V""^^: V^ = ifoff + ^io(L • B) + £:ii(Sp • 
B) + Ei2{Sd ■ B) + Eis{Se ■ B). In the adopted leading order approximation En = -4.2577 kHz/G, E12 = -0.6536 
kHz/G and E13 = 2.8025 MHz/G are expressed only in terms of the mass and magnetic moments of the particles. The 
values of Eiq were calculated using the variational non-relativistic wave functions of HD+ of Ref. [16!| , thus improving 
on the accuracy of the first such calculations [15]. 

The matrix of V^^ is evaluated in the basis set of vectors with definite values of the squared angular momenta 
F = Sp + Se, S = F + Sd, J = S + L and 2;-axis projection J^ of J. Except for the "stretched" states with the 
"extreme" values F = l, 5 = 2, J~L + 2 and Jz = ±J, Jz is the only exact quantum number; labeling the eigenstates 
of y* requires an additional index n. The corresponding eigenvalues AE'"^"'^^ represent the energy levels of HD+, 
defined relative to the "spinless" energies E^^, calculated as eigenvalues oi H^^ + V'^^^^. Since the hyperflne spectrum 
is only deformed but not rearranged by magnetic fields B below a few G, we take for n the set of quantum numbers 
{F, S, J) labelling the hyperfine states at B = 0. 

The hyperfine states of HD+ are split into sublevels distinguished with the quantum number J^ . The Zeeman shift 
may be approximated with the quadratic form 

The numerical values of ft,"^" and q"^"^^ calculated by the least square method and providing relative uncertainty 
below 10~^ for S < 1 G are given in TableU] Eq. ([T]) may be used to evaluate the Zeeman shift of transition frequencies 
and choose the less sensitive ones as candidates for precision spectroscopy. Stretched states are a special case: the 
Zeeman shift is strictly hnear, with /i''^" = (L Eio + E12 + (E'n + £'i3)/2)/(L + 2). 

Transitions. The probability per unit time for an electric dipole transition between a lower \i) = \vLnJz) and upper 
{f\ = {v'L'n'J'z\ states, stimulated by an oscillating electric field E(i) = E cos27ri/i of frequency ly, is (in units h = 1) 
7r2(5(i/-i//,)](/|VF]i)]2, where i^f^ = iyo + AE'''^'"''-^'- -AE''^"-'', i'q=E'"'^' -E"^ is referred to as "central frequency", 
M^ = — E.d is the interaction with the electric field, and d is the electric dipole moment operator. The matrix elements 
{f\W\i) are expressed in terms of the reduced matrix elements of d which, in the non-relativistic approximation, do 
not depend of the spin quantum numbers: dfi = dy'L' vL- Accurate numerical values of dyiL',vL have been calculated 
with the variational Coulomb wave functions of Ref. [la]; they agree with earlier results of Ref. [ijj at the 10"'^ level. 
Assuming statistical population of the hyperfine states, the observable spectrum of the transition (w, L) — )► {v' , L') can 
be put in the form ]E]^ J2fi'^i^ ~ '^/O^fi (^)j where 6 is the angle between B and E, the sum is over all pairs of 
states («, /) belonging to the hyperfine structure of the initial and final ro- vibrational states. The number of hyperfine 
lines may exceed 10'^, but most of them are weak. For magnetic fields B below a few G the spectrum is dominated 
by the "favored" components corresponding to transitions between states that satisfy AS* = AF — 0; the frequencies 
of the favored transitions lie in a band of width ^ 100 MHz around vq- Fig. [T] illustrates the Zeeman structure of an 
overtone ro-vibrational transition ((0, 1) — > (4, 2)), to be discussed further below. 

Of interest also are the radio-frequency (RF) magnetic dipole (Ml) transitions between spin sub-states of the 
same (vL) state, stimulated by an oscillating magnetic field B' cos 2T:i>t. The frequencies Vfi of the transitions with 
AF = ±1 are in the range 700 < Vfi < 1100 MHz, while z//i < 200 MHz for AF = 0. The exphcit expressions of the 
RF transition probabilities show that they are mainly due to the (B' • Sg) term in y™^s 

The two-photon spectrum may similarly be put in the form |E|^ ^^ ^ 5{v—Vfi/2) Tr^ [9). In transverse polarization 
{0 — 'k/2) each favored hyperfine component splits into 3 components with AJz — J'z — Jz = ~2,0,2 with a typical 
separation of the order of ^ 10^ kHz at i? = 0.5 G; each of these components acquires an additional super-fine 
structure with separations in the 10^ kHz range. In a longitudinal magnetic field {6 — 0) AJz — 0, and the lines are 
"super-fine" split only. 

Transitions between stretched states are allowed by selection rules, both as one-photon pure rotational {v = v') or 
ro-vibrational as well as two-photon transitions. 

Discussion. At low B, the Zeeman shift of the frequency of a single Zeeman component may be estimated using 
Eq. dl]): 5vfi{B) w {h' J'z — hJz)B + {q' — q)B^. The systematic uncertainty due to the external magnetic field may be 
minimized by selecting transitions with minimal Si'fi{B). We describe a few concrete transitions, for different cases 
of assumed experimental spectral resolution. 

1: Resolved hyperfine structure, unresolved Zeeman subcomponents. The level of spectroscopic resolution needed for 
improved measurements of mass ratios is about 10 kHz for an overtone vibrational transition and therefore requires 
experimental resolution of the hyperfine structure (with typical spacing of several MHz, see Fig.[T]), but not necessarily 
of the Zeeman sublevels. 

Among one-photon transition, a favorable case is the triplet of magnetic subcomponents \vL{FSJ)Jz) — 
01(010)0) — > ]42(011)J^). The J^ = — > J^ = subcomponent (parallel polarization) has a very small quadratic 
Zeeman shift (-18 Hz at 1 G), while the 7^=0 —?• J^ = ±l subcomponents (observable in perpendicular polarization) 
have dominant linear Zeeman shifts of 34.3 and —36.3 kHz, respectively, with an intensity-weighted shift of only 290 



Hz at 1 G (1.3 ■ 10~^^ in relative units). If the magnetic field direction is not optimal, the latter value sets the scale 
of the Zeeman systematic effect. A two-photon transition candidate is the ^-insensitive |00(011) J2) — > |20(011) J^) 
triplet (line M20 in Fig. [2]) with a weighed shift of —95 Hz at 1 G (—8 • 10^^'^). As fields below 1 G are feasible, such 
shift values will be well below the current and near-future theoretical uncertainties and thus will not be a limitation 
in comparisons of theoretical and experimental results. 

2: Resolved Zeeman subcomponents: the preferred experimental situation, in which the transition frequencies 
of individual subcomponents can be measured. Examples of one-photon transitions with very low magnetic field 
sensitivity include transitions with J'^ = Jz — with eliminated linear and suppressed quadratic dependence on 
B: \vL{FSJ)Jz) = |01(011)0) ^ |02(012)0} (pure rotational), |01(112)0) -^ |32(113)0), |01(010)0) -^ |42(011)0), 
101(012)0) -> |42(013)0), 102(012)0} -^ |43(013)0), with sensitivities 5v{lG)lvo « (3,0.1,-0.08,-0.2,-0.2) • lO^^^^ 
respectively. The sensitivity may also be weak in transitions with AJ^ ^ 0, e.g. in |02(011) — 1) -^ |43(012)— 2) with 
\6v{B)/vo\ < 8 • 10-" for B < 1 G. 

One-photon pure rotational or ro-vibrational transitions between stretched states occur as a doublet with purely 
linear magnetic shift: Svf^{B) = ±{L' E[q - LEw)B for J^ = ±{L + 2) ^ J^ = ±{L' + 2), where L,Eio {L',E[q) 
refer to the initial (final) state, and \L — L'\ — 1. This is about ±0.5 kHz for i? = 1 G. The mean of the frequencies 
is strictly independent of B in the adopted approximation. It should be possible to measure each doublet component 
independently and then compute the average, as in atomic optical clocks. 

Fig. [Dgives an example of a one-photon transition, (0, 1) — > (4, 2) at A « 1.4 /im. The hyperfine spectrum includes 
18 favored lines of which 10 are plotted. While most of them are quite sensitive to the external magnetic field 
strength, the Jz — ^ -^ J'z=^ subcomponents of M% and Mio are weakly sensitive to B. Line Mio also contains 
the moderately sensitive doublet mentioned above, and M4 contains the stretched states transitions observable at 
perpendicular polarization only. 

Two-photon transitions between stretched states are of metrological interest. In transitions with L' = L ±2 the 
doublet splitting is of order ±10Hz/G, since Eqi varies weakly with the vibrational quantum number, and the mean 
shift is zero. Again, it should be possible to measure each doublet component independently and then compute the 
mean. For two-photon transitions with L' = L — there is no splitting and no shift. Favourable two-photon transitions 
are (0, 0) — >■ (2, 0) at 2 x 5.362 /im and (0, 1) —^ (2, 1) at 2 x 5.368 /im, which have intermediate states whose detuning 
is not too large and therefore allow a reasonable excitation strength for realistic laser power from quantum cascade 
lasers. In addition, the (0,0) — ;■ (2,0) transition has a particularly simple spectrum, see Fig. [51 since the L — states 
have only four levels and the two-photon transition selection rules 18] allow only A J = AJz — transitions. This 
gives a relatively strong weight to the stretched states transition relative to the sum over all transitions (which are 
all equally strong) . Its frequency is given by a simple expression, see below. 

3. Comparison of ab-initio theory and experiment. When spectroscopy of HD+ is pursued with the goal of com- 
parison with QED calculations and for a determination of the particle mass ratio, it may be important to determine 
the central frequency vq since in the adopted approximation, unlike the spin corrections, vq includes the contribution 
of QED corrections of order 0{ma^) and 0(ma^ log a). We propose to combine the results of optical pectroscopy 
between vibrational levels {v, L) and («', L') with the results of RF spectroscopy within each of these levels to extract 
1^0- Indeed, the coefficients of the effective spin Hamiltonian Ei and the spin corrections to £'"^, can be expressed 
in terms of the frequencies of RF transitions within the hyperfine structure of the (v, L) state (and similar for the 
{v',L') state). Since the number Uh of linearly independent hyperfine transition frequencies exceeds the number iVe 
of non- vanishing Ei, the above relations can be resolved for Ei only if the radio- frequencies satisfy Nh — N^ com- 
patibility relations, that may serve to test the experimental accuracy. The value of vq is then expressed in terms of 
spectroscopic data and may be used to test the predictions of few-body bound state QED in next-to-leading orders 
in a. This approach is particularly simple in states with L = 0, since the only non- vanishing coefficients are E4 and 
E5. Denote the frequencies of the Jz—0 -^ J'z—O components at i? = of three RF transitions in the {v,L — 0) 
state as follows: i/i = ^'«o(iii)o,«o(ioo)o, ^^2 = i't,o(i22)o,«o(iii)o 1^3 = i^«o(ioo)o,«o(oii)o, and similarly, by i^[...,i^!^, the 
corresponding frequencies in the (v\0) state. Diagonalizing the hyperfine hamiltonian, we obtain £'4 ~ 1^2 + 1^3, 
£'5 = 2 (vi -\- 1^2)73, with compatibility relation 1/2(^^1 + ^^2 + I's) = 21^1 i^s (and similar for the (v' , 0) state). The spin 
energy of the stretched states is AE'"^'-^^^''^ = £'4/4 -I- E5/2 = i^i/3 + 7i^2/12 + 1^3/4. The transition frequency between 
stretched states is i^i,o(i22)2,ti'o(i22)2 = ^^o + Av, Av = 1/3 [v'^ — vi) + 7/12 (^ ^ ^2) + {vj, — i^3)/4. Thus, vq could be 
obtained by measuring one vibrational transition, and six RF transitions. The sensitivity to magnetic field comes only 
from Av, since the Zeeman shift of i^i,o(i22)2,D'o(i22)2 is strictly zero in the adopted approximation. For example, Au is 
shifted by only 55 Hz (5 • 10^^^) in the case w = 0, u' = 2 at _B = 1 G. This shift value could be taken as a conservative 
Zeeman uncertainty for vq, but it could be further reduced by correcting for the magnetic field in the trap; the field 
could be determined e.g. from a measurement of the Zeeman splitting of an appropriate magnetic-sensitive transition 
and using the theoretical magnetic sensitivities. 

^. Spectroscopic techniques. In order to achieve the spectroscopic resolutions discussed, a Doppler-free technique is 
required, and two methods appear suitable. Two-photon spectroscopy with counter-propagating beams of an ensemble 



of cooled, but not strongly confined molecular ions strongly suppresses the first-order Doppler broadening [19]. Second 
is quantum logic spectroscopy [20| on a single HD"*" ion. For both methods, techniques for strongly populating the 
lower spectroscopic state will be helpful |ll| . 

Summary. By evaluating the Zeeman effect for all experimentally relevant spectroscopies of HD+ we have shown that 
it will not be a limiting factor for the experimental accuracy if appropriate transitions and spectroscopic technique 
are selected. This is of particular relevance to the possible use of HD+ for setting limits to a hypothetical time- 
dependence of particle mass ratios, where a relative accuracy of 10^^^ is desired. For the determination of the ratio 
of electron to reduced nuclear mass by comparison of ab initio QED calculations and experimental frequencies, it is 
advantageous to determine the central transition frequencies vq since the uncertainties related to particle magnetic 
moments and nuclear electromagnetic structure will be significantly suppressed. We have shown how this can be 
achieved by combining two-photon and RF spectroscopy of states with low magnetic sensitivity, the latter feature not 
being affected by the higher-order relativistic mass and anisotropy corrections. 

Other systematic shifts, such as light, blackbody radiation and Stark shifts are currently being analyzed. 
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FIG. 1: Left: Strength Tr^ of ten favored hyperfine components Mi, . . . , Mio of the (0, 1) -^ (4, 2) one-photon transition line 
at 5 = (no ^-dependence) , labeled with the quantum numbers (FSJ) of the initial and final states. Right: Frequencies of the 
A Jz = components of Mi , . . . , Mio as function of the magnetic field B for 9 = 0. The Jz = — ^ J'z = components of Ms and 
Mio (short dotted lines) undergo a very small quadratic shift. The frequencies of the transitions between stretched states in 
transversal field (dotted lines) have a very small linear shift; their mean is independent of B. 
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TABLE I: Numerical values of the coefficients h"^" (upper line, in kHz/G) and q"'^"-'%J^ = (lower line, in kHz/G^) in 
the weighed least-square approximation (1) of the dependence of the Zeeman corrections to the hyperfine energy levels on 
the magnetic field B for the ro-vibrational states of the HD+ ion with L < 4 and v < A. The relative inaccuracy of the 
approximation does not exceed 5.10~ for B < 1 G, except for the F = 1 components of the states with L — \ where accuracy 
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